Self-trapping of strong electromagnetic beams in relativistic plasmas 
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Interaction of an intense electromagnetic (EM) beam with hot relativistic plasma is investigated. 
It is shown that the thermal pressure brings about a fundamental change in the dynamics - localized, 
high amplitude, EM field structures, not accessible to a cold ( but relativisic) plasma, can now be 
formed under well- defined conditions. Examples of the trapping of EM beams in self-guiding regimes 
to form stable 2D solitonic structures in a pure e-p plasma are worked out. 
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I. INTRODUCTION 

The problem of electromagnetic (EM) wave propaga- 
tion and related phenomena in relativistic plasmas has 
attracted considerable attention in the recent past. From 
the nonthermal emission of the high-energy radiation 
coming from a variety of compact astrophysical objects it 
has become possible to deduce the presence of a popula- 
tion of relativistic electrons in the plasma created in the 
dense radiation fields of those sources [1]. The principal 
components of these plasmas could be either relativistic 
electrons and nonrelativistic ions (protons), or relativis- 
tic electron-positron (e-p) pairs. 

Relativistic e-p dominated plasmas may be created in 
a variety of astrophysical situations. Electron-positron 
pair production cascades are believed to occur in pulsar 
magnetospheres [2]. The e-p plasmas are also likely to 
be found in the bipolar outflows (jets) in Active Galactic 
Nuclei (AGN) [3], and at the center of our own Galaxy 
[4]. In AGNs, the observations of superluminal motions 
are commonly attributed to the expansion of relativistic 
e-p beams in a pervading subrelativistic medium. This 
model implies copious pair production via 7 — 7 interac- 
tions creating an e-p atmosphere around the source. The 
actual production of e-p pairs due to photon-photon in- 
teractions occurs in the coronas of AGN accretion disks, 
which upscatter the soft photons emitted by the accre- 
tion disks by inverse Compton scattering. The presence 
of e-p plasma is also argued in the MeV epoch of the early 
Universe. In the standard cosmological model, tempera- 
tures in the MeV range (T - 10^"^ - IMeV) prevail up 
to times t = Is after the Big Bang [5]. In this epoch, the 
main constituent of the Universe is the relativistic e-p 
plasma in equilibrium with photons, neutrinos, antineu- 
trinos, and a minority population of heavier ions. 

Contemporary progress in the development of super 
strong laser pulses with intensities / ~ 10^^~^^W / cm? 
has also made it possible to create relativistic plasmas 
in the laboratory by a host of experimental techniques 



[6]. At the focus of an ultrastrong laser pulses, the elec- 
trons can acquire velocities close to the speed of light 
opening the possibility of simulating in the laboratory 
the conditions and phenomena that, generally, belong in 
the astrophysical realm [7] . 

Elucidation of the electromagnetic wave dynamics in a 
relativistic plasmas will, perhaps, be an essential tool for 
understanding the radiation properties of astrophysical 
objects as well as of the media exposed to the field of 
superstrong laser radiation. Although the study of wave 
propagation in relativistic plasmas has been in vogue for 
some time, it is only in the recent years that the nonlin- 
ear dynamics of EM radiation in e-p dominated plasmas 
[8] has come into focus. The enhanced interest stems 
from two facts : 1) e-p plasmas seem to be essential con- 
stituents of the universe, and 2) under certain conditions, 
even an ultrarelativistic electron-proton plasma can be- 
have akin to an e-p plasma [9] . 

Wave self-modulation and soliton-formation is, per- 
haps, one of the more interesting and significant signa- 
tures of the overall plasma dynamics. The existence of 
stable localized envelop solitons of EM radiation has been 
suggested as a potential mechanism for the production 
of micropulses in AGN and pulsars [10]. In the early 
universe localized solitons are strong candidates to ex- 
plain the observed inhomogeneities of the visible universe 
[11,12]. 

In e-p plasmas, there does not exist a general satisfac- 
tory theory for the generation of soliton like structures by 
ultrastrong high-frequency EM fields of arbitrary spatio- 
temporal shape. Relative complexity of the fully rela- 
tivistic equations (hydrodynamic or kinetic) has limited 
their solutions essentially to one-dimensional systems . A 
summary of the salient results is : in unmagnetized e-p 
plasmas, high-frequency pressure of the EM pulse pushes 
e-p pairs in the direction of its propagation thus creat- 
ing a density hump in the region of the field localization. 
The effective refractive index of medium decreases in this 
region and as a result localized soliton-formation is not 
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supported by the medium. (The increase in the refractive 
index due to the relativistic nonlinearity related with the 
particles' high-frequency motion is not enough to cope 
with the decrease caused by the above mentioned effect). 
In Refs.[ll, 13] it has been argued that locahzed solitons 
can be formed if the interaction between the EM field 
and acoustic phonons is taken into account- the envelope 
solitons propagating with subsonic velocities may, then, 
emerge. In magnetized e-p plasma, larger classes of soli- 
ton solutions are possible-typical examples may be found 
in Refs.[14]. However, it is conceivable that soliton solu- 
tions obtained in a one dimensional formulation will turn 
out to be unstable in higher dimensions. 

It is, therefore, a matter of utmost priority that we 
explore the possibility of finding stable multidimensional 
soliton solutions in e-p plasmas. Dynamics of 3D enve- 
lope solitons of arbitrary strong EM fields in e-p plasma 
with a small fraction of heavy ions has been analyzed in 
Ref.[12]. It was shown that, in a transparent e-p plasma, 
EM pulses with L\\ « L± (where L|| and L± respec- 
tively,the characteristic longitudinal and transverse spa- 
tial dimensions of the field), may propagate as stable, 
nondiffracting and nondispersing objects (light bullets) 
with large density bunching. It was further shown in Ref. 
[15] that these bullets are exceptionally robust: they can 
emerge from a large variety of initial field distributions, 
and are remarkably stable. In these studies, the EM field 
is pulse-like with longitudinal localization much stronger 
than the transverse; the localization is brought about by 
the charge separation electric field (usually absent in a 
pure e-p plasma) created by the presence of a small pop- 
ulation of ions. 

In the present paper we explore another mechanism 
for localization - we will show that the in the pure e- 
p plasma, the thermal pressure can provide the confin- 
ing "glue" just as the charge separation electric field did 
for an e-p plasma contaminated with ions. We will also 
deal with a complimentary manifestation of the radia- 
tion field- the "beam" {L± « L||) rather than the pulse 
{L± >> X||). Assuming the plasma to be transparent to 
the beam, and applying a fully relativistic hydrodynam- 
ical model, we demonstrate the possibility of beam self- 
trapping leading to the formation of stable 2D solitonic 
structures. The high-frequency pressure force of the EM 
field ( tending to completely expel the pairs radially from 
the region of localization) is overwhelmed by the thermal 
pressure force which opposes the radial expansion of the 
plasma creating conditions for the formation of the sta- 
tionary self-guiding regime of beam propagation. 



II. BASIC EQUATIONS 

In this paper the word "relativistic" connotes two dis- 
tinct regimes: the plasma becomes relativistic when ei- 
ther the directed fluid velocity approaches the speed of 
light or the thermal energy per particle is of the order 
of, or larger than the rest mass energy. Since both these 



paths to relativity are encountered in the astrophysical 
as well as laboratory plasmas ( produced and accelerated 
by intense laser pulses), we will investigate a fully rela- 
tivistic hydrodynamical model. 

If the velocity distribution of the particles of species 
a (= e,p,i, ... where e,p, and i denote respectively elec- 
trons, positrons and heavy ions (protons)) is taken to 
be a local relativistic Maxwellian, the hydrodynamics of 
such fluids is described by [16]: 

^{W^UakWa) - ^Pa = -F'^'Jak (1) 
OXk OXi c 

where = [70, 7q:Uq,/c] is the hydrodynamic four- 
velocity with Ua as the three- velocity , 7^ = (1 — 
u\/(?)~^/'^ is the relativistic factor, Jak is the four- 
current, F**^ is the electromagnetic field tensor, and Wa 
is the enthalpy per unit volume 

Wo. = "^mocC^Go. f (2) 

7a V^OaC / 

Here tooq and are the particle invariant rest mass and 
temperature, respectively, n„ is the density in the labo- 
ratory framc^ of the fluid of species a. The pressure Pa = 
UaTaha, and Ga{za) = (l/^a)/ A'2 (l/z« ) , where K2 
and K3 are, respectively, modified Bessel functions of the 
second and third order and = Ta/moaC^ . The factor 
Ga{za) has the following asymptotes, Ga ~ 1-|-5zq/2 for 
Za «l { non-relativistic) and Ga ~ 4^^ for >> 1 ( 
highly relativistic). 

The set of equations (l)-(2) may be written in the stan- 
dard form : 

%(moaC^G„7a) - — = eaUct ■ E (3) 
at Ua ot 

%(GaPa) + — VP„ = e„E + ^(u„ X B) (4) 
at Ua c 

where = "/ainoaUa is the hydrodynamical momen- 
tum, E and B are the electric and magnetic fields, and 
da/dt = d/dt + Ua ■ V is the comoving derivative. The 
hydrodynamical velocity Uq, and the relativistic "fa are re- 
lated to the momentum by the standard relations: = 
Pa/moa7a and 7„ = (l-|-p„^/mQ„c^)^/^. It is interesting 
to note from Eqs. (3)- (4) that the fluid inertia is modified 
by the temperature; the expression = moaGa{za) 

denotes the effective mass of the particle. For ultrarela- 
tivistic temperatures (Tq, >> moaC^), the effective mass 
turns out to be M^fj = ^Ta/c^ » nioa- Thus the par- 
ticles "forget" their rest mass and the plasma tiuns into 
a kind of "photon" gas. If an ultrarelativistic plasma is 
in thermodynamical equilibrium with the high-frequency 
photon gas {huj ^ T), one should also take into account 
the radiation pressure Pji = aT^ {a = w/45h^c^) [17]. 
In this paper this effect will be neglected. We must also 
bear in mind that for extremely relativistic regimes, the 
model Eqs.(3)-(4) fail to adequately describe the plasma 
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dynamics since heavy particle production has been ne- 
glected. This shortcoming will impose an upper limit on 
the temperature for the validity of the model.. Note that 
in the context of early universe, the epoch in which the 
e-p plasma is dominant, has a characteristic temperature 
Ta ^ IMeV and Meff ~ 4moe. Since the particle masses 
are just a few times larger than their rest mass at these 
temperatures, the e-p plasma can still be considered as a 
two component fluid rather than a photon gas. 

The equation of state directly follows from the self- 
consistency of Eqs.(3) and (4). Taking the scalar product 
of Eq.(4) with and comparing it with Eq.(3), we can 
derive dalnPa/dt = ZadaGa/dt. After straightforward 
manipulation, the equation can be easily integrated to 
yield 



Pa = C exp(2;aG„) 



(5) 



where the arbitrary constant C must be defined by the 
initial state. Using = naTa/'fa Eq.(5) reduces to the 
adiabatic equation of "state" : 



la K2{Za) 



exp(— = const 



(6) 



For nonrelativistic temperatures, Eq.(6) yields the usual 
result for a monoatomic ideal gas {rf^/Ta^'^ = const, 
where = na/la is the density in the rest frame 
of fluid element) and for ultrarelativistic temperatures 
one obtains the equation of state for the photon gas 
(n^/T^ = const). We would like to emphasize that the 
7a factor appearing in Eq.(6) is related to the coherent 
or directed motion of fluid elements whose origin may 
lie either in an initial macroscopic flow or in the motion 
imparted by intense EM radiation. 

To complete the picture we must couple the plasma 
equations of motion with Maxwell equations: 

cV X B = — 47rJ (7) 



cVxE = -- (8) 



where 



V • E = 47rp (9) 



V-B = (10) 



p = ^eocna, J = ^eanaUa (11) 



are, respectively, the charge and current densities. The 
system of Eqs.(3)-(11) along with the continuity equation 
( for each species) 



dna 



V • (riaUa) = 0, 



(12) 



represents a closed set of equations which describe prop- 
agation of EM radiation in relativistic multicomponent 
plasmas. 

The above system can be manipulated further to reveal 
interesting structural properties. To begin with, Eqs.(5)- 
(6) can be cast in the form 



mpaC 
Ja 



(13) 



which, when substituted into Eq.(4), converts it to 

d 

g^{GaPa) + Woac2V(G„7a) = e^E + [u„ X (14) 



where 



J2« = ^B + Vx (G,pJ 
c 



(15) 



is the so called generalized vorticity . Taking the curl of 
Eq.(14), we find that the evolution equation for 



V X [Ua X ila] 



(16) 



is of the standard vortex dynamics form. Although the 
system of Eqs.(14)-(16) can be traced to early publica- 
tions (see for instance Ref.[18]), their consequences are 
yet to be fully worked out. An immediate consequence, 
for instance, is the appropriate equivalent of Kelvin's 
theorem- the flux of generalized vorticity is frozen- 
in through a comoving area. 

The system yields the following set of relativistic 
Beltrami- Bernoulli equations for equilibrium [19]: 



Gala + 



const 



(17) 



(18) 



where are constants and (j) is the scalar potential. The 
relevance of these equilibria for astrophysics is the subject 
matter of a forthcoming paper. 

For the current effort, we apply Eqs.(14)-(16) for wave 
processes in an unmagnetized plasma. From Eq.(16) it 
follows that if the generalized vorticity is initially zero 
(fla = 0) everywhere in space, it remains zero for all 
subsequent times. We assume that before the EM radi- 
ation is "switched on", the generalized vorticity of the 
system is zero. In this case the equation of motion may 
be written as: 



d_ 

dt 



Ha + moaC Vr„ = e^E 



(19) 



where the temperature dependent momentum Ua and 
r„ are defined by: 



Ua — GaPo 



Ta = Gala = V Gl + {Ua/niOaCy 



(20) 



(21) 
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The condition of vanishing generalized vorticity connects 
Ua with the magnetic field: 



B = V X Ua 



(22) 



It is remarkable that in Eq.(19) the magnetic part of the 

Lorentz force is formally absent; this fact greatly sim- 
plifies analytical manipulations. It is equally remarkable 
that our equations which describe the dynamics of a hot 
relativistic plasma are structurally similar to equations 
used in the theoretical treatment of different aspects of 
ultrastrong laser interaction with a cold plasma [20] . This 
similarity becomes even more evident when we study the 
interaction of short EM pulse with relativistic electron- 
ion plasmas. If the pulse is assumed to be shorter than 
the characteristic time for ion response (i.e. inverse of 
ion Langmuir frequency), the ion motion may be ignored, 
and the electric field may be found from the electron part 
of Eq.(19), 



eE = - moeC^VFe 

ot 



(23) 



Substituting this expression into Poisson's Eq.(3) (which 
now reads as V • E = 47re(noi — rig), where riQi is the 
equilibrium ion density) we find the electron density 



rie 



1 



1 



d 



(24) 



where = (47re^noi/moe)"^/^ is the plasma frequency. 
Using Eqs.(22)-(24), Eq.(7) reduces to : 



c^V X V X He + 



5^ 



+ moeC 



dt 



-+ 



2 He 



1 



1 



9. 



c2 

^^^^ 



= 



(25) 



moewf dt 

which, along with the equation of state {zg = moeC^/Te) 

^a^e Zq 



K2{Ze) 



exp(— Ge-2e) = const, 



(26) 



constitutes the simplified system to which the entire set of 
Maxwell and relativistic hydrodynamic equations for an 
the electron -ion plasma have been reduced. An equation 
similar to Eq.(25) was derived in the cold plasma limit 
in Ref. [21]. However, there arc important differences : a) 
Due to the temperature dependence of Gg in Eq.(21), the 
factor Fe and the momentum He are no more related by 
simple relations as they do for a "cold" case, and b) to 
incorporate the temperature variation in the system we 
must add the equation of state (26). 

Though Eqs.(25)-(26) form a more complicated system 
than what we have for the cold plasma, we believe that 
many results obtained in the cold plasma limit can find 
appropriate analogies in the hot relativistic- temperature 
case. Detailed studies of the nonlinear dynamics of the 



electron-ion plasma is beyond the intended scope of the 
current paper and some new results will be presented 
separately elsewhere. In the next part of the current 
paper, we concentrate on the nonlinear dynamics of EM 
beams in e-p dominated plasmas. 



III. THE ELECTRON-POSITRON DOMINATED 
PLASMA 

In this section wc apply our general formulation to 
the problem of self-trapping of EM beams in pure e-p 
plasmas with relativistic temperatures. For notational 
convenience, wc replace the subscripts (c) and (p) by 
superscripts (— ), and (+). Wc assume that the equilib- 
rium state of the plasma is characterized by an overall 
charge neutrality n~ = n+ = rioo, where and n+ 
are the unperturbed number densities of the electrons 
and positrons in the far region of the EM beam localiza- 
tion. In most mechanisms for creating e-p plasmas, the 
pairs appear simultaneously and due to the symmetry of 
the problem it is natural to assume that = T+, = Too, 
where and T+ are the respective equilibrium temper- 
atures . 

We shall assume that for the radiation field of in- 
terest, the plasma is underdense and transparent, i.e., 
e = We/w <C 1, where ui is the mean frequency of EM 
radiation and u>e = {^ire^noo/mQeY^'^ is the plasma fre- 
quency. Since both species are mobile, the c-p dynamics 
can not be reduced to just one vector equation similar 
to Eq.(25). We will display the entire set in terms of 
potentials ( the Coulomb gauge V • A = will be used). 



E = - 



IdA 

c~dt 



■V</), 



B = V X A, 



(27) 



and the dimensionlcss quantities t — ujt, f = (a;/c)r, 
= T^/moeC^, A = eA/(moec2), (j) = e^/moeC^, 
= Ti^ /{moec), and = n^/uao- Suppressing the 
label- "tilde", we may arrive at the dimensionless equa- 
tions. 



dt 



(28) 



^-^^+Ft^'^- ''^^^ - J") = (29) 



A(j) = e^(n~ - n+) 
V- A = 



dt 



(30) 
(31) 

(32) 
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=n±/r± and r± = s/iG^Y + (n±)2. The To = const. Note that from Eqs.(30) and (33), we can 



with J* = 
species equation of state is: 



r±/(T=' 



where 



K2il/T^)T^ 



exp[G=^/T±] 



(33) 



(34) 



Of various techniques that could be invoked to investigate 
Eqs.(28)-(34) to study the self-trapping of high-frequency 
EM radiation propagating along the z- axis, we choose 
the method presented in the excellent paper by Sun et al. 
[22]. The method is based on the multiple scale expansion 
of the equations in the small parameter e. Assuming that 
all variations are slow compared to the variation in ^ = 
z — at, we expand all quantities Q = (A, </>, 11^, 71=*=, ...) 
as 



Q = Qoi^,xi,yi,Z2) +(-Qi{^,xi,yi,Z2) 



(35) 



where {xi,yi, Z2) ~ {ex,ey,e^z) denote the directions of 
slow change, and ai = {a^ — ~ 1. We further 

assume that the high-frequency EM field is circularly po- 
larized, 

Ao_L = ^ (x -h iy)^ exp(i^/ a) + c.c. (36) 

Here A is the slowly varying envelope of the EM beam, 
X and y denote unit vectors, and c.c. is the complex con- 
jugate. We now give a short summary of the steps in the 
standard multiple- scale methodology (Ref. [22]). To the 
lowest order in e, we obtain the following. The transverse 
(to the direction of EM wave propagation z) component 
of Eq.(28) reduces to 



ta, 



0_L 



(37) 



and for the longitudinal components we get: 



an, 



02 



Equations (29)-(31) yield d^V±4>o = 9|0o = ^^^^o = 
0, where V_l is the perpendicular Laplacian in {xi,yi). 
These relations imply that and Aqz do not depend on 
the fast variable ^. For the self-trapping problem, we 
can assume that Aq^ = IIq^ = [22]. From Eq.(38), and 
from the lowest order continuity Eq. (32), we obtain : 
(J^Fq = d^riQ = 0, i.e., Fq and Uq also do not depend 
on the fast variable ^. 

To the next order (in e), the transverse component of 
Eq.(28) reads : 

- ^"'4 + V^F± = T(-a)^ T V^<^o (39) 



Averaging over the fast variable ^ we obtain Vj_rg = 
TV_l(/'o yielding the trivial solution (pQ = and Fq = 



deduce that 



no 



and Tq = Tq 



Thus, as one would expect, the low frequency motion 

of the e-p plasma is driven by the ponderomotive pres- 
sure (~ IIoj.^) of the high-frequency EM field and this 
force, being same for the electrons and positrons, does 
not cause charge separation. It is also evident that due 
to the symmetry between the electron and positron flu- 
ids, their temperatures, being initially equal, will remain 
equal during the evolution of the system. The relation 
between the EM field and the temperature can be found 
using the equation Fq = const obtained above. Using 
Eqs.(36)-(37) and by choosing the const by requiring that 
at infinity A — > and Tq Too, we derive 



G^iTo) = G^iT^)-\A\' 



(40) 



It follows from Eq.(40) that the present hydrodynami- 
cal model, which describes the nonlinear waves in e-p 
plasma, is valid for |Ap/G^ < 1. In our opinion the ori- 
gin of this restriction lies in the inadequacy of the basic 
model, and is not solely due to a failure of the perturba- 
tion technique used above. When this condition is vio- 
lated, the EM waves are overturned and they will cause 
multistream motion of the plasma (i.e. wave breaking 
takes place). In such a situation, one must resort to ki- 
netic description for studying the nonlinear wave motion. 
Notice however that the function G{Tao) ^ 1 if Too 
but rapidly increases with increase of Too thus providing 
room for l^jmax to reach from weak to relativistic values. 

We are now ready to deal with the equation for the 
slowly varying envelope A of the EM beam. To the lowest 
order in e, one finds from Eq.(29) 



oi- 



32 A, 



0_L 



-ViAo±-2 



d^dz2 



o± ,2!^Ao^ = (41) 



Go 



F or deriving thi s equation, we used the relation Fq 
^/G^iT^fTW = Goo, and 



no (To) = 



f{To) 

/(Too) 



(42) 



which follows from Eq.(33). Substituting Eq.(36) into 
Eq.(41) we find: 



oz Go 



■[l-n(To)]A = 



(43) 



where subscripts for variables {z2,xi,yi) are dropped for 
simplicity. We also assumed without loss of generality 
that (a^ — l)/e'^a'^ = 2/Goo, which in dimensional units 
coincides with the linear dispersion relation of the EM 
wave in an e-p plasma, namely: ui^ = liJ^jG^ -\- fc^c^ 
provided that a = uj/kc. 

Thus, the dynamics of EM beams in hot relativistic 
e-p plasma has become accessible within the context of a 
generalized nonlinear Schrodinger equation (NSE) (43). 
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IV. THE SELF-TRAPPED BEAMS IN E-P 
PLASMA 

In this section wc seek the localized 2D soliton solu- 
tions of Eq.(43), and analyze the stability of such solu- 
tions. Making the self-evid ent re-n ormalization of vari- 
ables z zG oo , r± — > r_L-\/Goo/2, Eq.(43) can be writ- 
ten as: 



i— + V]_A + ^A = 
oz 



(44) 



where ^ = 1 — no (To) represents the generalized non- 
linearity. The companion equation (40) can be viewed 
as a transcendental algebraic relation between Tq and 
1^1 p. Thus we conclude that ^l* is a function of \A\'^ 
(^f = ^-d^p). We note that Eq.(44) can be written 
in the Hamiltonian form iA^ = 5H/SA*, where H = 
/rfr^[|Vi^|2 - F(|A|2)] and F{t) = Jl;^it)dt. The 
Hamiltonian structure implies that Eq.(44) conserves the 
Hamiltonian H in addition to the power ("photon num- 
ber") N = J dr^\A\^. 

Unfortunately, it is not possible, in general, to derive 
an explicit analytic relation ^ = ^{\A\^) for arbitrary 
value of Too- Some qualitative deductions readily follow. 
Equation (40) shows that the presence of EM radiation 
reduces the temperature Tq. Since df{To)/dTQ > 0, from 
Eq.(42) we conclude that the plasma density is also re- 
duced in the region of the EM field localization which is 
in accordance with adiabatic motion of the plasma. For 
higher strength of the EM field, a complete expulsion of 
plasma i.e. plasma cavitation can take place (no ^ 0); 
this has been predicted in Ref . [23] . Thus the nonlinearity 
function "if shows a saturating character with the increase 
of EM field strength (note that present model is valid pro- 
vided |^|2/(G^ - 1) < !)• To illustrate, we exhibit in 
Fig.l a plot of \1/ versus |Ap for Too = 0.1. One can see 
that the nonlinearity function indeed saturates at high 
intensity. For small temperatures, we can even obtain 
an analytic expression for the function ^P. Remember- 
ing To < Too, assuming Too ^ 1, and by using Eq.(42) 
along with the asymptotic expansions Go(~ 1 + 5To/2) 
and /(« Tq), we derive for the nonlinearity function: 



1 - 



\Al 
5To, 



3/2 



(45) 



Equations (44)- (45) admit a stationary, nondiffracting 
axially symmetric solution of the form A/^/8T^ = 
U{r) QyiY>{iXz) where r = [x^ +1]^)^^'^ and A is the nonlin- 
ear wave-vector shift. The radially dependent envelope 
U (r) obeys an ordinary nonlinear differential equation: 

d?'U IdU T/rr9Nrr „ /.^N 

—^ + -—-XU + ^{U^)U = Q (46 
dr-^ r dr 

where 'J' = 1 — (1 — C/^)^/^. This equation corresponds to 

a boundary value problem with the boundary conditions: 
U has its maximum Um at r = 0, and [/ ^ as r ^ oo. 




FIG. 1: The nonlinearity function ^ versus \A\^ for Too = 0.1. 



We remind the reader that it has been shown in a semi- 
nal paper of Vakhitov and Kolokolov [24] that such solu- 
tions exist for arbitrary saturating nonlinearity functions 
^, provided that the eigenvalue A satisfies < A < "im, 
where '^m is a maximal value of the nonlinearity function. 
Equation (46) admits an infinity of discrete bound states 
characterized by j = 0, 1,2... zeros at finite r. We con- 
sider only the lowest-order nodeless solution of Eq.(46), i. 
e. "ground state" that is positive and monotonically de- 
creasing with increasing r. In the asymptotic region the 
solution must decay as Ur ^oo ^ exp(— ■\/Ar)/v^. Our 
nonlinearity function ^ has a maximum = 1 found 
at Urn{^ l)i i- ^- at the maximally allowed strength of 
the field. As a consequence the upper bound of the prop- 
agation constant Ac must satisfy Ac < ^m- Numerical 
simulations show that the amplitude of the ground state 
solution Um = U{r = 0, A) is a growing function of A 
(see Fig. 2) and it acquires its maximum value (= 1) at 
A = Ac « 0.29. The solution represents a trapped, local- 
ized EM solitary beam. The beam becomes wider at low 
amplitudes. 

The stability of the solitonic solutions can be investi- 
gated using the criterion of Vakhitov and Kolokolov [24]- 
the soliton is stable against small, arbitrary perturba- 
tions if dN/dX > 0, where N is the power of the trapped 
mode: 



iV(A) = 27r / dr rU^{r,X) 



(47) 



In Fig. 3 we plot the numerically obtained solutions of N 
for various A. Since the curve has positive slope every- 
where, the corresponding ground state solution is sta- 
ble for < A < Ac. Notice that the power of the 
solitary beam always exceeds a certain critical value 
N > Nc ^ 7 .8 . We also know that N must be bounded 
from above (iV < iV„ « 10.5). 

For arbitrary temperatures, explicit form of = 
^'(|A|^) can not be found. However, due to its satu- 
rating character, solutions with properties similar to the 
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FIG. 2: Nonlinear dispersion relation: the amplitude Um as 
a function of A, the eigenvalue. The plasma temperature 
Toe « 1- 
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FIG. 3: The beam power N versus A (Too << 1). 



small temperature case (which can be explicitly solved) 
could be expected. Using relations (34), (40) and (42), 
we numerically find a stationary solution of Eq.(44) for 
arbitrary T^o ■ For convenience we use following represen- 
tation of vector potential A/A^ = U{r) exp{iXz), where 
Ac = (G^ — 1)^/^. Though the maximum value of U 
is still restricted by the condition < Um < 1, the am- 
plitude of vector potential A^ can reach a considerable 
value. For ultrarelativistic temperatures, Too >> 1, we 
have Ac = y/TET^ » 1 and since < Am < Ac the soli- 
ton solution with ultrarelativistic strength of EM field 
is possible. Here we present results of simulations for 
Too = 1 (i.e. Too ~ O.bMeV ). In Fig.4 wo plot the 
amplitude of the ground state solution Um versus the 
propagation constant A. The solution exists provided 
< A < Ac(« 0.22). The profiles of the field U{r) the 
plasma density no(r) and the temperature To(r) are ex- 
hibited in Fig. 5 for A = 0.1. One can see that in the 
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FIG. 5: Normalized EM field Z7, plasma temperature To and 
density no versus r for Too = 1- 

region of field localization, the plasma temperature and 
density is reduced. Similar plots could bo obtained for 
all allowed values of A. When A Ac, plasma cavitation 
takes place, i.e. at r = the plasma density and temper- 
ature tends down to zero. Appearance of zero temper- 
ature is not surprising since the corresponding region is 
the "plasma vacuum"; all particles are gone away. 

The dependence of iV on A is presented in Fig. 6. One 
can see that the curve N = N{\) has a positive slope 
and according to Vakhitov and Kolokolov criterion, the 
corresponding solitary solutions are stable against small 
perturbations. 

The detailed dynamics of arbitrary field distribution 
must be studied by direct simulations of Eq.(44). We 
can learn much , however, from the recent extensive elu- 
cidations of the dynamical properties of the solutions of 
NSE with saturating nonlinearity. It seems that the gen- 
eral features of evolution are not sensitive to the details of 
the saturating nonlinearity (see for instance [15, 25] and 
references therein). For all such systems the beam will 
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FIG. 6: The beam power versus A (Too = 1). 

monotonically diffract if tlie beam power is below a criti- 
cal value {N < Nc), and it will be trapped if N > Nc and 
the Hamiltonian H < is negative. In the latter case, 
the beam parameters will oscillate near the equilibrium, 
ground state values. This oscillations monotonically de- 
crease with increase z due to the appearance of the radi- 
ation spectrum. For larger z, the oscillations are damped 
out, and the formation of the soliton in its ground state 
takes place. If the initial profile of the beam is close to the 
equilibrium one, the beam quickly reaches the ground- 
state equilibrium, and propagates for a long distance 
without distortion of its shape. The initial beam, even 
when its parameters (i.e. amplitude, effective width and 
phase) are quite far from equilibrium, will either focus or 
defocus to the ground state, exhibiting damped oscilla- 
tions around it. As a consequence the beam reaches an 
equilibrium with its final power slightly smaller than the 
initial. Such an evolutionary scenario may not hold for 
very intense beams with N >> Nc] the beam may then 
break up into filaments due to a modulation instability. 
However, each filament, will tend to evolve towards its 
own equilibrium state corresponding to the power it car- 
ries. Thus, the ground-state equilibrium seems to be an 
attractor. 

Our own studies indicate that Eq. (44) , with the nonlin- 
earity particular to the problem at hand, reproduces the 
general expected behavior described above. However we 
find that the soliton formation requires the initial beam 



power to be in the range Nc < N < Nm- For N > Nm, 
the multistream motion of the plasma prevents the sys- 
tem from settling in a steady state. 



V. CONCLUSIONS 

We have investigated the nonlinear propagation of 
strong EM radiation in a relativistic, unmagnetized two- 
fiuid plasma. The treatment is fully relativistic- in the 
coherent or directed motion as well as in the random or 
thermal motion of the plasma particles. The assumption 
that prior to the switching of the field-plasma interaction, 
the generalized vorticity is zero, greatly simplifies the sys- 
tem of relativistic fluid equations. In particular, in the 
electron-ion dominated plasma, under well defined con- 
ditions the system of Maxwell -fluid equation ( Eqs.(25)- 
(26)) turns out to be structurally similar to the one ob- 
tained for a cold plasma. Consequently we would expect 
that results already established in cold plasma limit can 
find appropriate analogy in the hot plasma case. 

We presented a somewhat detailed study of EM beam 
propagation in transparent e-p plasmas. Applying a re- 
ductive perturbation technique, the system of relativistic 
Maxwell- fluid equations is reduced to a 2D nonlinear 
Schrodinger equation with a saturating nonlinearity. We 
found that if the strength of the EM fleld amplitude is 
below the wave breaking limit, the beam can enter the 
self-trapped regime resulting in the formation of stable, 
self-guided 2D solitonic structures. The beam-trapping 
owes its origin to the thermal pressure (which opposes 
the ponderomotive pressure) - Naturally such trapping 
regimes are not accessible in the relativistic but cold 
plasma limit. In the region of beam trapping, the plasma 
density as well as its temperature is reduced and under 
certain conditions these parameters can be reduced con- 
siderably (i.e. plasma cavitation takes place). 

The fact that relativistically hot e-p plasmas are ca- 
pable of sustaining high amplitude localized structures 
of high amplitude electromagnetic fields should be a re- 
sult of considerable importance to an understanding of 
the complex radiative properties of different astrophysi- 
cal objects where such plasmas are known to exist. 
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